Spin-Pseudospin Coherence and CP 3 Skyrmions 
in Bilayer Quantum Hall Ferromagnets 
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We analyze bilayer quantum Hall ferromagnets, whose 
underlying symmetry group is SU(4). Spin-pseudospin co- 
herence develops spontaneously when the total electron den- 
sity is low enough. Quasiparticles are CP 3 skyrmions. One 
skyrmion induces charge modulations on both of the two lay- 
ers. At the filling factors = 2/m one elementary excitation 
consists of a pair of skyrmions and its charge is 2e/m. Recent 
experimental data due to Sawada et al. [Phys. Rev. Lett. 
80, 4534 (1998)] support this conclusion. 

The quantum Hall (QH) effect is a remarkable macro- 
scopic quantum phenomenon in the two-dimensional 
electron system jjj. Attention has recently been paid 
to quantum coherence in QH systems. The kinetic and 
Coulomb Hamiltonians have the spin SU(2) symmetry. 
Its spontaneous breakdown leads to a spin coherence and 
turns the system into a QH ferromagnet. The effective 
Hamiltonian is the SU(2) nonlinear sigma (NLer) model 
P). Quasiparticles are CP 1 skyrmions ^-^|. 

In this paper we analyze skyrmion excitations in bi- 
layer QH (BLQH) ferromagnets. The lowest Landau 
level (LLL) contains four energy levels corresponding to 
the spin and layer (pseudospin) degrees of freedom. The 
SU(4) symmetry underlies the BLQH system. Its spon- 
taneous breakdown leads to a spin-pseudospin coherence. 
The effective Hamiltonian is the SU(4) NLcr model 0. 
Quasiparticles are CP 3 skyrmions Q. They have two 
characteristic features: (A) One skyrmion induces charge 
modulations on both of the two layers. The main part 
of the activation energy is the capacitive charging en- 
ergy. (B) One elementary excitation consists of a pair of 
skyrmions at v = 2/m with m odd. It carries the charge 
2e/m. Our theoretical analysis accounts for recent ex- 
perimental data due to Sawada et al. |(],|7j . Throughout 
the paper we use the natural units h = c = 1 . 

QH ferromagnets: We analyze skyrmion excitations at 
the filling factor v = 2jrpo / eB± = 1 /m with m odd. We 
use an improved composite-boson (CB) theory B, which 
is proposed based on a suggestion due to Girvin et al. 
P|. The advantage of the scheme is a direct connection 
between the semiclassical property of an excitation and 
its microscopic wave function. We start with a review of 
monolayer QH ferromagnets ||] . The analysis of BLQH 
ferromagnets is its straightforward generalization with a 
replacement of SU(2) by SU(4). 

The kinetic Hamiltonian for planar electrons with 
mass M in the perpendicular magnetic field B± is 



where P 9 
A 



-idj + eAj is the covariant momentum with 
j - 7j£j k x k B±; £i 2 = -e 2 i = 1 and en = e 2 2 = 0. 
The electron field ^(a;) is a two-component spinor made 
of the spin-up (ip^) and spin-down (ip±) field. The sym- 
metry group of this Hamiltonian is U(2)=U(1)<B>SU(2). 

When the Zeeman energy is small, a spin coherence 
develops spontaneously. This is described by introducing 
the two-component CB field by the formula 



$(as) = e 



(2) 



where A(x) is the auxiliary field obeying 'V 2 A(x) — 
2-Km[p(x) — po\; p(x) = W(x)9(x) is the electron 
density. The phase field 9 (a;) attaches m units of 
flux quantum 27r/e to each electron via the relation, 
£ijdidjQ(x) = (2w/e)mp(x). The effective magnetic 
field is 

B cB (x) =B±- EijdidjQix) =B X - (2ir/e)mp(x). (3) 

It vanishes, (# e ff)g = 0, on the ground state at v = 1/m. 
Substituting (||) into (Q), the kinetic Hamiltonian reads 



11 h = 2M 1 d2x ^^ V - 



iVy)(V X +iVyMx), (4) 



where Vj = —idj — (ejk + iSjk.)dkA(x) is the covariant 
momentum. We have defined &(x) = $^ (x)e 2A ^ x \ with 
which p(x) — ^(x)^>(x) — &(x)Q(x). An analysis of 
the Lagrangian shows that the canonical conjugate of 
<p a (x) is not iif a '(x) but i<p a ^(x). 
We decompose the CB field as 



&(x)=e- A M<l>(x)n{x), 



(5) 

\J p{x) and the 



with the U(l) component cj)(x) — e ix ^ x '^ 
SU(2) component n(x): It is the CP 1 field subject to 
the constraint, n'(x)n(x) = 1. The spin density is 

S a (x) = ^p(x)s a {x), s a (x) = n){x)\ a n{x), (6) 

where A a are the Pauli matrices. 

At sufficiently low temperature it is reasonable to focus 
our attention to physics taking place within the LLL. 
The Hilbert space Hlll is constructed by imposing the 
LLL condition so that the kinetic term (^) is quenched. 
It has a simple expression in terms of the CB field, 



(v x + iV y )Hx)\e) = - J-_$(x)|s 



0. 



(7) 



The complex number is z = (x + iy)/2lB with £b the 
magnetic length. Hence, the wave function for composite 
bosons is analytic and symmetric in all N coordinates, 



1 



n[z] = (o\${x 1 )---$(x N )\&). 



(8) 



The wave function for electrons is &[x] — fl[z]&i,j<i[x], 
where 6ln[^] is the familiar Laughlin wave function. 
Here, [x] = (xi, x 2 , ■ • • , x N ) and [z] = {z\, z 2 , ■ ■ ■ , z N ). 
Any excitation confined to the LLL is described by a 
choice of the analytic spinor factor Q[z]. 

We analyze the CB theory semiclassically, where the 
bosonic field operator is approximated by a c-number 
field. It follows from m) that the wave function is 



e[x]= n<$(a; r ))6 LN [x] 



(9) 



where ($>(x)) is analytic. If the Zeeman energy is ne- 
glected, the ground state is determined by minimizing 
the Coulomb energy. The one-point function ($(x)) is a 
constant vector pointing to an arbitrary direction in the 
SU(2) space, as implies a spontaneous breakdown of the 
SU(2) symmetry. Actually, a small Zeeman interaction 
fixes this direction so that (s a (x)) — S az , or 



(¥> T (»)g = VPo, (<^O))g = 0. 



(10) 



The ground-state wave function is given by (^) with ([lO|) . 
The semiclassical LLL condition follows from (||), 



(<p a {x)) = e- A{x) e lx{x] ^{x)n a {x) = LU a (z) 



(11) 



where various fields are classical ones. This is solved by 
X(x) = and n a (x) = uj a (z) / VEJ^"! 2 '■ The soliton 



— V 2 \np(x)-p(x) 



equation |10| follows trivially from (|ll|), 

Po = vQ(x), 
where Q(x) is the topological charge density, 
g( ; r) = i-V 2 ln(^K(z)| 2 ). 

The lightest skyrmion on the ground state ([Toj ) is 

(fHx)) sk y = (p l (x)) sky = Ky^, 



(12) 



(13) 



(14) 



with which the wave function is given by ( |9|). For a large 
skyrmion (n 1), the soliton equation (|12|) is solved 
iteratively and agrees with the familiar expression S, 



g(x) = p(x) - po - vQ(x) = - 



7r [r 2 -I- 4(k£ b ) 



212 ' 



(15) 



The skyrmion is reduced to the vortex in the limit k — > 0. 

The excitation energy of one skyrmion consists of the 
exchange energy E ex , the electrostatic energy Ec and the 
Zeeman energy Ez ■ Minimizing their sum, we determine 
the skyrmion size k, the skyrmion energy -E s k y and the 
flipped-spin number A^ spin as follows 



-1/3 



E, 



ky- 



^spin— 2i/k 2 In 



2n 



(16) 
(17) 
(18) 



Here, p s — ^e 2 /(16v27re^ B ), Eq — e /e£b and g — 
g* pbB I Eq, The parameter f3 represents the strength of 
the Coulomb energy; it is calculated as (3 — 37r 2 /64 for 
a sufficiently large skyrmion in an ideal planar system. 
However, an actual skyrmion is not large and there will 
also be a correction from a finiteness of the layer width. 
We use (3 as a phenomenological parameter. We fix it as 
(3 = 0.24 by requiring the skyrmion spin -/V sp i n to agree 
with the experimental data due to Barrett et al. U at 
v = 1: ~ 3.7 at B = 7.05 Tesla, where K ~ 1.0 and 

spin ~ 5.3 (k ~ 1.10) at B = 3 



spin 

0.015. Note that N s 



Tesla and A^ spin ~ 2.7 (k ~ 0.96) at B = 15 Tcsla. Our 
results are consistent with the previous ones ||. 

The excitation energy of a skyrmion-antiskyrmion pair 
will be given by A = 2E sky — r o g- set with a sample de- 
pendent offset r o g- set . In Fig.^we have fitted the experi- 
mental data due to Schmeller et al. M based on formula 
(|l7| ) with (3 = 0.24, where an appropriate offset r o ff set 
is used for each curve. The theoretical curves reproduce 
all data remarkably well. 
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Zeeman energy g*fJB tot (e£) 
FIG. 1. Theoretical curves versus experimental data due 
to Schmeller et al.[5] for the activation energy in three sam- 
ples SI1, QW1 and QW2. There are two curves for one sam- 
ple (QW1) but with different mobilities. The offset r o ff se t 
increases as the mobility decreases. 

BLQH systems: We proceed to analyze BLQH states 
at v = 1/m and 2/m with m odd. There are three ex- 
perimental techniques to elucidate various states. The 
first one is to change the total electron density po- 
By increasing po the interlayer separation d is effec- 
tively increased compared with the magnetic length, 
d/ts oc t/pq. Hence, as po — > oo the BLQH system 
at v = 2/m will be decomposed into two independent 
v = 1/m monolayer systems; as po — > an interlayer 
coherence will develop spontaneously, which has been 
argued |11 lQ,h3| at v = 1/m and will be argued at 



2 



v = 2/m in this paper. Consequently, we expect a phase 
transition at v = 2 /m between these two phases but not 
at v = 1/m. These two phases are clearly distinguish- 
able by using the second technique, i.e. by applying gate 
bias voltage. We can control the density difference ao be- 
tween the two quantum wells, where ao = (Po — Poll Po 
with Pq the density in the layer a. Only the coherent 
state is stable [|ll] against an arbitrary change of cr . All 
these features have been experimentally confirmed in re- 
cent works due to Sawada et al. at v = 1 and 2. The 
coherent state at v = 2/3 has not been observed by them 
presumably because of a poor sample quality. The third 
technique is to tilt the sample with the perpendicular 
magnetic field B± fixed. In the high-density data jjj the 
activation energy is found to increase at v = 2 and 2/3 
as normally as in the monolayer system: Indeed, we can 
fit the data by the monolayer skyrmion formula (|l7|). 
We conclude that elementary excitations are monolayer 
skyrmions. In the low-density data @] it is found to 
decrease anomalously at v = 1 and v = 2, as is the phe- 
nomenon discovered by Murphy et al. at v = 1: It 
is a behavior intrinsic to the interlayer coherent state in 
the BLQH system. 

Spin-pseudospin coherence: We analyze elementary 
excitations in the coherent state of the BLQH system. 
The electron field ^(x) has four components ip 1 ^ , i/) 1 *, 
ip 2 ^ and ip 2 ^- , where the layer is indexed by 1 and 2. 
The kinetic Hamiltonian is given by (|l|), whose symme- 
try group is U(4)=U(l)<g>SU(4). When the interlayer and 
intralayer Coulomb energies are nearly equal and domi- 
nate the system, we expect a spin-pseudospin coherence 
to emerge. 

Such a new phase is described in terms of the CB field 
defined by (|^). The CB field is decomposed into the U(l) 
and SU(4) components by (||). Here, n(x) is the CP 3 
field. The group SU(4) is generated by the hermitian, 
traceless, 4x4 matrices A a , a = l,2,---,15, normalized 
as Tr(A a A & ) = 26 ab . They are the generalization of the 
Pauli matrices in case of SU(2). The SU(4) spin density 
is given by (||) with such A a . 

The Coulomb energy is decomposed into two terms, 

E C = IJ d 2 xd 2 yV+(x - y)g±(x)g ± (y), (19) 

where V ± (x) = (e 2 /2e)(\x\~ 1 ±(\x\ 2 + cP)- 1 ^) with the 
interlayer separation d; g + = g(x) = p(x) — p ; Q-{x) = 
p 1 (x) — p 2 (x) -/?o + /°o with p 1 (x) — p^(x) +p 1J -(cc) and 
p 2 (x) = p 2 ^(x) + p 2 ^(x). The Coulomb energy Eq, pos- 
sessing the SU(4) symmetry, is the driving force to realize 
the QH system. The term Eq describes the capacitive 
charging energy between the two layers. It vanishes in 
the limit d — > 0. 

The Hilbert space Hlll is defined by the LLL condi- 
tion (0). In the semiclassical approximation the wave 
function is given by © at v < 1. The SU(4) spm- 
pseudospin coherence is shown to develop spontaneously, 



precisely as the SU(2) spin coherence is. The ground 
state at v < 1 is given by [jl0| 

<^) g = VaTVTT^, (^) g = o, 

(^ 2 T) g = VP^VT^, (^ 2i )g = 0, (20) 

where p\ = |po(l + c ) and pi = |po(l - <?o). It is 
convenient to use a new set of CB fields, 

</ T = V(l + ^o)/2^ 1T + V(l-^o)/2^ 2T , 

^= V(l - a )/2/T _ ^(1 + ^/2^, (21) 

and a similar set for the spin-down fields. They are re- 
duced to the symmetric and antisymmetric fields at the 
balanced point (ao = 0). We call them the "bond" and 
"antibond" fields. The ground-state value (^0|) is trans- 
formed into ($) g = v / po(l, 0, 0, 0) in terms of the new 
fields, i> = (tp s 1 ,tp SL ,ip A ^ ,ip Al ). 
The tunneling energy is 

E T = -±A SAS ^l-a 2 J d 2 x[g s _(x) - g A {x)], (22) 

where g^_(x) = p s (x) — (p s (x)) g with p s = <p s ^(p s ^ + 
(p s lt(p s -l and similar equations for g_(x). The tunneling 
gap is (1 — a 2 ) 1 / 2 Asas on the state (|20|). 

A skyrmion excitation flips in general spins and pseu- 
dospins. It is a CP 3 skyrmion described by 

(^)sky = y/po(z,Ki,K2,Ka), (23) 

with constant parameters «j. Its classical configuration 
is determined by (0) - (|l|), and @ with k 2 = J2 a K l- 
For definiteness we assume hereafter that the tunneling 
energy is much larger than the Zeeman energy. (For in- 
stance, Asas/<?*M-B^ — 4 at B = 5 Tesla in the sample 
of Ref. ||.) Then, we have k = Ki ^ 0, «2 = «3 = 0. It 
is identical to the CP 1 skyrmion ([l4]) in the spin space. 
At the balanced point the skyrmion size, the skyrmion 
energy and the flipped-spin number are given by the 
same formulas as ( |T6| ) ~ (^|), where (3 now depends on 
the layer separation d. Here, we concentrate our atten- 
tion to its dependence on the imbalance parameter ao. 
The term involving ao is only the charging term ( |l9| ) for 
v < 1. We give a numerical estimation of the activation 
energy at v = 1 in Fig.|| by using the sample parame- 
ters (d = 23lA, £ B = 120A) of Ref. §. The theoretical 
curve explains the experimental data quite well with a 
reasonable skyrmion size k ~ 1. 

BLQH ferromagnet at v — 2: A caution is needed to 
analyze the BLQH system at v = 2 since one Landau 
state contains two electrons. We attach one unit of flux 
to each electron and transform it into the CB field by 
formula (^|) with m = 1. The effective magnetic field is 
not given by (||) but by 

2B cS = 2B ± - EijdidjQix) = 2B ± - {2n/e)mp(x). 

(24) 
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It vanishes on the homogeneous ground state at v = 
2/m. Due to the fermi statistics the wave function is the 
antisymmetric product of two wave functions (0), 

6[X] = Y[[{§l(x r )) (8 (®2(x r )) - (* 2 (av)> ® ($l(x r ))] 
r 

xSlnW 2 , (25) 
where ($>\(x)) and ($2(2;)) are analytic and satisfy 

X>?(*)>| 2 =X>?(*)>| a , (26) 

a a 

as follows from the semiclassical LLL condition flll] 
10 
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FIG. 2. Theoretical curves versus experimental data due 
to Sawada et al.[6]. The skyrmion charge is e at v = 1 and 
2e at ^ = 2. The discrepancy of the data for large |<to | may 
indicate that genuine CP 3 skyrmions are excited there since 
the tunneling gap (1 — oo) 1 / 2 Asas becomes smaller. The 
dotted curve is for a would-be skyrmion carrying e at v — 2. 

When Asas 3> <?*Mb-^j the spin-up and spin-down 
"bond" states are filled. Hence, the ground state is given 
by (E5h with a set of two constant CB fields, 



-(1,0,0,0), 



($5 



-(0,1,0,0), (27) 



in terms of the "bond" and "antibond" fields. This 
might be identified with the canted state Q for Asas 3> 
g*lisB. A skyrmion excitation flips pseudospins, or in- 
duces a coherent tunneling excitation. It is described by 
(pB|) with a set of two CB fields, 



(3>2)sky= \fpo/2(K[,Z, k' 2 , k' 3 ), 



with K 2 = J2 a K a = J2 a K 'a- H consists of two CP 3 
skyrmions (pf), and the skyrmion charge is 2e. We em- 
phasize that there exists no skyrmion composed of one 
CP 3 skyrmion at v = 2 because of the constraint (26). 

An estimation of the excitation energy of the skyrmion 
(p8|) is straightforward. We concentrate our attention 
to its dependence on the imbalance parameter <7o. The 
terms involving <jq are the charging energy (|l9| ) and the 
tunneling energy (S3). The charging energy increases 



while the tunneling energy decreases as oo increases. We 
give a numerical estimation in Fig.^ by using the sample 
parameters (d = 23lA, l B = 228A, A S as =6.8K) of 
Ref. ||. The vortex limit (k ~ 0) gives a best fitting of 
the experimental data because of a large tunneling gap 
(Asas =6.8K). We have also given a theoretical curve 
for a would-be skyrmion carrying charge e at v — 2 by 
using the same parameters, where the charging energy 
and the tunneling energy are found to cancel each other 
almost completely (Fig.|^). 

The driving force of the SU(4) spin-pseudospin coher- 
ence is the Coulomb exchange energy arising from the 
SU(4)-invariant Coulomb term Eq in (|l9|). Provided the 
exchange energy is dominant, it is obvious that the SU(4) 
coherence develops also at v — 2/m with the skyrmion 
charge 2e/m, and at v = 6, 10, 14, • • • with charge 2e. 
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